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SUl^RY 

We  first  solve  the  equation  dX  +  aXdt  =  dN,  where  dN 
represents  a  Poisson  process,  and  then  generalize  to  a  Levy 
process.  Finally,  we  solve  a  linear  partial  differential  equa¬ 
tion  DX  =  dL  in  strong  distribution,  meaning  that  the  second 
member  dL  is  a  distribution  process,  generalization  of  Levy 
process  on  F.  The  results  are  then  applied  to  wave  propa¬ 
gation  in  underwater  acoustics,  and  spatial' correlation  is 


determined. 


INTRODUCTION 


The  theory  of  linear  random  fields  is  well  developed;  for  example  S. 
Badrikian  and  S.  Chevet  [2]  consider  the  cylindrical  measures  and  the  associ¬ 
ated  linear  random  field.  Some  applications  to  economy  and  industry  are 
presented  in  S.  Ben  Soussan  [3] .  We  consider  the  linear  random  process  X 
which  is  the  solution  of  a  linear  partial  differential  equation  formally 
written  DX  =  dT,  where  dT  is  a  random  process;  in  the  simplest  case  dT  is 
white  noise.  The  equation  Dp  =  0  in  the  last  chapter  of  this  paper  will  be 
the  propagation  equation  of  the  pressure  p  in  deep  sea  water.  We  consider 
only  linearized  propagation  equation  from  general  ones  (cf.  Poirde  [11]);  our 
equation  is  thus  an  approximation  in  the  sea  medium.  In  D.  de  Brucq  [4],  the 
second  member  of  the  equation  is  a  Wiener  measure  or  Gaussian  measure,  denoted 

4 

by  dW,  defined  on  S(F  )  the  space  of  indefinitely  differentiable  functions  on 

4 

F  decreasingly  quickly.  The  second  member  describes  the  random  approximation. 
This  equation  is  written 


Dp  =  Ap 


1  ^  ^  b  3 

■  2  ^^2  2  at 

c  3t  c 


Ap  =  dW 


where  A  is  the  Laplacian  in  F  .  The  solution  in  the  sense  of  strong  distri- 

4 

bution  meaning  satisfies  for  any  f  of  S{F  )  : 

Tf.?' 


p(f)  -  /  T 


dW 


where  F  is  the  Fourier  transform  on  F 


and  T  its  inverse  and  where  A  is  the  function  defined  by  the  equality  F*D  = 

2  2 

AF.  In  this  particular  case  A(u),x)  =  l2(u){x  -[k(u))-iY(u))]  }  where  k  is  the 

4  * 

wave  number,  y  is  the  absorption.  By  hypothesis,  for  any  f  and  g  of  S(F  )  ,  the 

2  "  2 

correlation  of  the  noise  dW  is  r(f,g)  =  E(W(f)WCg))  =  a  JfgdX  with  o  a  normali- 

4 

zation  constant  and  dA  Lebesgue  measure  on  F  .  The  spatial  correlation  s  at 
a  given  frequency  v  =  ^  of  the  pressure  p  is  given  (cf.  Th  III-3-7)  by; 
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s  (u), !,) 


2  2 


IT  a 


X 


'YH  .  .  . 

e  sin  k)t 

—  'X  — ; - 

Y  k£ 


here  8,  is  the  distance  between  the  points  where  the  pressure  is  measured;  by 
physical  approximation  may  be  taken  equal  to  1. 

For  1=0,  the  experimental  spectral  density  function  of  the  random 
fluctuations  of  the  pressure  in  deep  sea  satisfy  to  a  good  approximation  the 
relation 


s (u,0) 


2  2 


TT  a 


I 

Y 


In  the  ocean,  point  processes  and  Poisson  processes  dN  appear  to  be  more 
accurate  than  Gaussian  processes  to  describe  the  random  sources  of  the  noise. 
We  generalize  as  much  as  possible  to  introduce  a  Levy  linear  process  dL  as 
second  member  of  the  equation.  The  expression  s(u),Jl)  does  not  change  with 
that  extension. 

In  the  first  section,  we  explain  the  method  with  the  linear  differential 
equation 

dX  +  aXdt  =  dN  a  >  0  (1) 


where  dN  is  the  centered  Poisson  process.  We  compute  the  characteristic  func¬ 
tion  of  the  X  process  so  that  the  Gaussian  or  the  Poisson  laws  of  the  forcing 
term  can  be  separated.  Not  withstanding  the  correlation  function  r  of  X  and 
in  the  stationary  situation  the  spectral  density  power  function  are  given. 

In  the  second  section,  we  introduce  a  Levy  process  L  with  stationary  incre¬ 
ments  and  recall  a  decomposition  theorem  for  this  indefinitely  divisible  pro¬ 
cess.  Then  we  solve  equation  (1)  using  L  as  forcing  term.  We  are  now  in  the 
situation  to  consider  spatio-temporal  problems.  In  the  third  section,  we  use 

4 

general  theorems  (A.  Badrikian  [1])  to  construct  a  measure  L  on  S(F  )  general¬ 
ization  of  the  Levy  processes.  We  solve  in  the  sense  of  strong  distribution 
the  equation  DX  =  dL.  The  case  where  L  is  Gaussian  measure  is  known  (D.  de 
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Brucq  and  C.  Olivier  fS]).  Poisson  measure  is  well  known.  (J.  Neveu  flO])  but 
DX  =  dL  with  L  Poisson  measure  is  new,  a  fortiori  with  L  our  generalization  of 
a  Levy  process. 

Finally,  we  apply  the  results  to  the  propagation  equation  and  we  obtain  the 
spatial  correlation  at  a  given  frequency. 

Some  classical  notations  and  results  will  be  useful.  We  will  write  the 
Fourier  transform 

F(f)(z)  =  J  dz  where  the  function  f  is  integrable 

r" 

and  where  <z,Z>  is  the  scalar  product  in  .  In  the  spatio-temporal  appli¬ 
cation  of  propagation,  for  physical  reasons 

F(f)  (oj,5,ri,C)  =  /  (t,x,y,z)dtdxdydz 

R^ 

Fourier  transformations  are  isomorphism  of  distribution  spaces  S(r”)  ,  SCR^^)  . 

The  distributions  space  SCR*')  is  nuclear  and  countably  semi-normed; 
these  notions  are  for  example  defined  in  I.M.  Gelfand  and  N.  Yu.  Vilenkin  [6]. 

We  denote  by  S  this  space  and  S'  its  topological  dual,  the  tempered  distri¬ 
bution  space.  The  topology  of  S  is  defined  by  the  semi-norms: 

Vk,  NcR,  s  (f)  ^  sup  l(l+|zr)Vf(z)| 

z  r" 

!a|<N 

we  have  the  topological  inclusion  P  c  S  c  with  V  the  space  of  indefinitely 
derivable  functions  with  compact  support  in  R**  . 

If  r  is  a  real  or  complex  function  on  r"  ,  continuous  in  zero  and  of 
positive  definite  type  then  the  Bochner  theorem  asserts  the  existence  of  a 
bounded  positive  measure  p  such  as 

Vt  e  r"  r(t)  =  /e^''^dp(v) 

The  random  processes  of  the  first  two  paragraphs  are  defined  on  a  probabil  ity 
space  (fi,A,P).  We  denote  by  R"  the  tribe  of  the  borelians  of  r"  . 
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I.  SOLUTION  OF  THE  EQUATION  dX  +  aXdt  =  dN 

We  take  a  to  be  strictly  positive  real  number  and  N  to  be  a  Poisson 
process. 

I . 1  Poisson  Processes 

Consider  a  family  of  random  variables  N{Tj,T2)  with  ^^>^2  e  F  ’ 

Each  has  Poisson  law  and  represents  the  number  of  discontinuities  of  a  random 
phenomenon  in  the  interval  [tj,T2]-  a(T)  be  a  non-decreasing  function 

defined  by  the  relation 

E(W(Tj^,T2))  =  /  d0(T)  . 

"l 

We  suppose  a  continuous  on  F  so  a  is  almost  bounded.  There  exists  a  sequence 
(T^.neF)  of  intervals  such  as 

T  and  U„T„  =  F 

n  n+1  n  n 

b)  J  do(2)  <  ® 

T 

n 

The  Poisson  process  N  of  parameter  a,  is  the  fcimily  of  random  variables: 

N(t)  =  N(0,t)  for  T  >  0  and 

N(t)  =-N(0,t)  for  T  <  0  . 

The  centered  Poisson  process  N  is 

Vt e  F  N(t3  =  M(t)  -  a(T)  and  we  write 
dN  =  dW  -  do 

The  solution  X  of  the  equation 

dX  +  a  X  dt  =  dN  is 
Vt  e  F  X(t)  =  /^„e'^^^"'^^dNCT) 

The  integral  is  the  almost  sure  limit  of  e"®^^~^^dN(T)  when  t  converges 
to  . 
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We  suppose  that  ^^do(T)  <  <»  and  also  that  2a(t  ^ 

Since  X  is  a  second  order  process,  straightforward  computation  gives 

y 

ft  -2a(t-T).  ,  .  I  ft  -a(t-T),  .  ,1 
E(X  (t))  =  J_^e  Mo{t)  +  J_^e  ■'do(T) 

and  the  variance  is 

E(X^t))-[E(X(t))]^  =  /^^e'^^^’''^^do(T). 

For  the  centered  Poisson  process  N,  the  solution 

Vt  e  F  X(t)  =  ^^^~^^dN(T)  is  centered  and 

E(X(t)^)  =  /^^e'^^f'''''^do(T) 

if  N  has  stationary  increments  then  a  is  a  Haar  measure  on  Rand  do  =  X^dx 

2 

with  dx  hebesgue  measure  and  constant.  For  X^  =  1,  E(N(t)  )  =  |t|  and 

E(X(t)^)  =  If  we  change  the  Poisson  process  N  to  a  Gaussian 

2 

process  W  such  as  E(W(t))  =  0  and  E(W(t)  )  =  t  then  the  solution  is  centered 
and  has  the  same  variance. 

The  characteristic  functions  of  N(x)  and  of  N(x)  are: 

'^W(x)(u^  =  E(e^“^*'’^b  =  expIo(x)(exp  iu-1)] 

'*'N(x)(u^  =  E(e^^’*-'^b  =  exp[o{x)  (exp(iu)- 1-iu)] 

1 . 2  Probability  Law  of  the  Solution  X 

Lemma  1.2.1:  The  complex  function  i|»  defined  by 

VfeL^(R  ,R,da)  CE)  =  /  [e^^^^^-l-i  f(z)]do(z)  is  continuous  on  L^(R  ,R,da) 

^  R 

2 

Proof:  For  every  f  of  L  (R  ,R,do),  we  have 

kj(f)  I  i  i/f^(T)do(x) . 

2 

Then  if  (f^)  converges  to  zero  in  L  (R  ,R,da)  then  converges  to  zero 


in  IE  . 


□ 
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Theorem  1.2.2:  There  exists  a  linear  centered  process 

N  =  (fi,A,P,(N(f))  ,  ,F  ,R) 

l/(R  .R.da) 

2 

with  characteristic  function  Vf  e  L  (F  ,R,do) 

A  iN(f). 

4>^  =  E(e 

Proof:  The  application. 

1.  L^(F.R,dT)  ^  N(t  )  -  N(t  )  £  L^(n,A,P) 

1^1 ^  ^ 

is  isometric  and  conserves  the  scalar  product.  The  family  1,  ,T,,T-e 

ITi,T2J  1  i- 

is  total.  We  note  N  the  unique  extension.  So  for  every  f  in 
2 

L  (F  ,R,da),  N(f)  is  a  centered  second  order  random  variable.  We  have  to 
compute  its  characteristic  function.  For  any  expression 

f  =  a^  1|^^  ^  1  in  L^(F  ,R,do) ,  we  have  N(f)  =  [N -N(t^ ) ] 


and 


'N(f) 


i  a  (N(t.  ,)  -  N(t  )) 

=  Ilog  E(e  :>  J  ) 

j 

ia . 

=  Ife  ^-1-i  a^]  [a(T^^p-a(Tj)] 

ila.l.  . 

=  /  {e^  ^'-^^-l-i  f(T)}da(T)  . 

F 

2 

By  continuity,  the  expression  is  valid  for  any  f  of  L  (F  ,R,da) 


□ 


Corollary  1.2.3;  The  characteristic  function  (f)^  of  the  process  X  solution  of 
the  equation 


dX  +  aXdt  =  dN 


satisfies  Vu,  t  e  F  . 


=  exp 


^  /!oo[®^P  e“^^^‘'^^-l-iu  e"®*^^''^^]da(T) 


Proof:  The  expression  of  X(t)  is 


with 

in  L^(R  ,R,da)  then  uX(t)  =  N(u  1  (•)  e“^*^^“^) 

[-“.t] 

Corollary  1.2.4:  The  covariance  function  r  of  the  process  X  solution  of  the 


equation 


dX  +  aXdt  =  dN  satisfies  Vt^.t^f^  F 


r(tj.t2)  =  E(X(t^)X(t2))  =  /_  e 


da(T) 


Proof:  As  XCt)  =  N(1  (•)  e  we  have 

[-".t] 

-a(t.-T)  -a(t»-T) 

rct-.t-)  =  /  1  (t)  e  1  (t)  e  da(T)  P 

^  F  [-».tj]  [-.t2] 

Corollary  1.2.5:  In  the  stationary  case,  the  covariance  function  r  is  equal  to 

^^^1*^2^  =  ^  e  ,  tj^,t2  in  F  ,  and  the  power  spectral  density  function 

is 

’'W  =S  *  ■ 

a  +v 

Proof:  In  that  case  da  =  X^dr  and  r(tj,t2)  “  ^  ® 

Fourier  transform  of  which  is 


I  /•+“  iVT  0  -a  x  ,  1  0 

TT— J  e  T5-  e  '  'dx  =  •= - = — ^ 

2Tr-'-“  2a  2-11  2  2 


II.  SOLUTION  OF  THE  EQUATION  dX  +  aXdt  =  dL 


We  generalize  the  forcing  term  to  be  a  Levy  process.  We  have  to  give 
some  essential  properties  of  these  processes  (T.  Hida  [9]). 

II . 1  Classical  Properties  of  the  Levy  Processes 

Definition  II.  1.1:  A  process 

L  =  (12, A, P,  (L(t))j^  ,  F  ,R)  is  a  levy  process  if  L  satisfies 
the  properties 

(a)  L  has  independent  increments. 

(b)  L  is  continuous  in  probability. 

(c)  L  has  trajectories  almost  surely  continuous  to  the  right  an^  .,ich  limit 
to  the  left. 

We  will  show  how  to  obtain  the  generality  of  such  processes. 

Proposition  II. 1.2:  Let  N  =  {Wjil=[a,bj  c  F}  a  family  of  processes.  For 
each  interval  1,  the  process  is  Poisson  such  that 

(a)  Vt  £  F  EfN^(t)]  =  t  n(I)  where  n  is  a  positive  measure  on  F  satisfying 
,2 

/  dn(z)  < 

|z|>0  1+z 

(bl  For  lj»l2  disjoin  (Ij^  1^  =  0)  then  the  Poisson  processes  Wj  are 

independent . 

(c)  For  every  partition  (I.  ;k  e  12 )  of  I=[a,b]  then  Vt  c  F  N  (t)  =  ^  N  (t) 

^  ^  kfF  ^k 

almost  surely. 

As  I  ->■  Wj(t)  is  a  measure  with  value  in  L°(fi,A,P),  we  use  the  notation 
.  The  process 

Vt  £  F  L^(t)  =  lim  /  ^  [z  dn(z)] 
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is  then  defined  and  is  a  Levy  process  with  stationary  increments.  The  measure 
n  is  the  Levy  measure  of  the  process  L^ .  The  convergence  with  p  is  almost 
surely  uniform  in  t. 

Remark  11.1.3:  Every  Poisson  process  is  a  Levy  process  and  the  vector  span  of 

such  processes  is  composed  also  of  Levy  processes.  We  need  only  to  be  able  to 

pass  through  the  limits.  Let  e  IN  }  a  sequence  of  independent  processes 

with  (Vt  e  R)  ^  o,  |(t)|  <  “;  is  'l  z  W  a  Levy  process?  If 
k  e  hi  keU 

(Vte  R)  I  z^a.(t)  <  “  then  ^  exists  and  is  a  Levy  process. 

k«N  ^  ^  keiN  ^  ^ 

Another  construction  is  possible.  Let  W  be  a  Poisson  process  and  (Hj^.k  ZL  ) 
a  sequence  of  independent  real  random  variables,  independent  of  W.  All  the 
discontinuities,  the  jumps  of  N,  are  equal  to  +1,  we  change  the  value  +1  of  the 
k^^  jump  in  +nj^,  random  in  R  and  we  obtain  a  Levy  process  (T.  Hida  [9]1.  We 
recall  the  P.  Levy's  decomposition  theorem: 

Theorem  1 1. 1.4:  Let  L  be  a  Levy-process  with  stationary  increments.  Then 
there  exist  two  constants  a,  6  of  R  ,  a  Wiener  process  W,  and  a  Levy  process  1.^ 
such  that  Vt  ^  R  L(t)  =  a  +  6  W(t)  +  L^Ct).  The  decomposition  is  unique. 

The  forcing  term  of  the  equation  dX  +  aXdt  -  dL  will  be  a  Levy  process  so 
that  we  generalize  the  Gaussian  and  the  Poisson  case. 

1 1 . 2  Solution  of  the  Equation  DX  *  aXdt  =  dL 

We  assume  the  Levy  process  L  with  stationary  increments  and  the  Levy 
2 

measure  n  satisfying  J  z  dn(z)  <  We  can  suppose  the  process  L  centered. 

R 

Proposition  1 1. 2.1:  The  process 

(Vt  e  R)[,^(t)  =  lim  J  it^^dz*-^^  '  ]  • 

p-H»  p>|z|>  - 

2 

where  W  is  the  Poisson  process  of  Proposition  II. 1.2  and  J  z  dn(z) 
uz 


<  OO 


is  a  centered  Levy  process  with  stationary  increments. 


Proof:  The  limit  J  i ^dM^  exists  ff.  Ilida  f9],  p.  39)  and  J 


f  ^cdn(z)  +  f^zdn(z)  +  f  ,zdn(z)  the  first  two  terms  converge  noting  for 

P  ^  \>\z\>^ 

'  ‘  P 

1  12  2 

example  that  [ /'  2dn(z)  j  <.  /_  z'^dn(z)  ^  /  z  dn(z)  <  ”  .  For  every  p  t  F, 

-P  -P 

f;[/  zW  (t)]  =  t  /  ,z  dn(z)  so  passing  through  the  limit 

P'>hl>^  P>h^l>i 

h(hjt))  =  0. 

2 

Lemma  11.2.2:  The  complex  function  define  by  Vf  e  L  (P  .R.dx) 

1)-Cf)  =  /  f  fexp  i  z  f(T)-l-i  z  f(T)]dn(z)dT 

RlzhO 

is  continuous  on  L^(F,R,dT). 

2 

Proof:  For  every  f  of  L  (F  ,R,dT),  we  have 

U,(f)|<.J  ilzl^lf(T)  |"dn(z)dT 

F  |z|>0 

<  /  |f(T) l^dT  /  lz|^dn(z)  n 

F  lz|>0 

2 

Theorem  1 1.2. 3:  If  /z  dn(z)  <  ■=»,  there  exists  a  linear  centered  process 


=  (n,A,P,(L^(f))  2 


,F  ,R) 


l.(F,R,dT) 

2 

with  characteristic  function  Vf  e  L  (F  .R.dr) 

i  L  (f) 

=  E(e  )  =  exp  i|;2(f) 


Proof:  We  can  suppose  J  z^dn(z)  =  1.  For  f  =  1.  ,,  the  characteristic 

F  ^  ’  J 


function  of 
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L^(t)  =  lim  J  j  (t) -tzdn(z) ]  is  CVu  f  P) 

p-v+co  p>|2|>i  ’ 

j''  J  [exp  iuz  1  (t1  -  -  iuz  1  (t)  dn(z)dT  . 

l^ho  10. t]  [O.t] 


As  J_  J  |i  u  z  1  (t)  |‘^dn(z)dT  <  “>  ,  the  random  variable  L  (t)  is  second 
■”  |z|>0  [O.t] 


order  and  E(L^(t)  )  =  t 


z  dn(z)  =  t.  The  same  is  true  for  t  <  0  and 


z  >0 


f  =  1 .  The  process  L^(t)  is  a  Levy  process;  the  map 


'[ti.Tjl  '  I-  -  ‘■c'V-'-cl’l)  ' 


(a.A.p) 


is  isometric  and  conserves  the  scalar  product.  The  family  ^  e  P 

is  total.  We  note  the  unique  extension.  So  for  every  f  in 
2 

L  (P  .R.dr)  L^ft)  is  a  centered  second  order  random  variable.  To  obtain  the 
characteristic  function  of  L^(t).  we  use  a  plain  extension  of  the  proof  of 
Theorem  1.2.2.  H 


Proposition  1 1. 2. 4:  Let  Vt  R  L  (t)  =  lim  J  j  [z  dn(z)  and 

^  p->co  p>|z|>l 

2 

siippose  /  z  dn(z)  <  «>  then  the  process  X  solution  of  the  equation 
lzl>0 

dX  +  a  Xdt  =  dL^  with  a  >  0  is  centered,  second  order,  stationary  with  covariance 


-alt^-tj 


Vt 


t^c  R  rctj.t^)  =  E(X(t2)XCtj))  = 


2a 


z  dn(z] 


z  >0 


Proof:  The  solution  X  satisfies 

Vt  (T  P(t)  =  /^^e"^^^"^^dL^(T) .  But  with  a  >  0  the  function 

^  is  element  of  L^(R  .R.dr)  then  X(t)  =  L^(l|^__^  . 

We  need  only  the  covariance  function  for  L  : 

c 
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)L.flt  ^  ,c 


f  r  I  -a(tj-T)  -a(t^-T)  2 

i  J  ,  t  1®  ^1-00  t  1®  “  dn(z)dT 

K  U  h-O  ^  ’  1  '  '  ’^2' 


2a  J  I  z I >0 


z^dnCz) 


Corollary  11.2,5:  The  power  spectral  density  function  of  the  solution  X  is 
Vv  t  R 

=  2V-2-^  /,  ,  ^^<ln(z) 

a  +v  1 2  I >0 

The  second  order  properties  are  the  same  for  every  Poisson,  Gaussian,  Levy 

centered  processes  and  for  the  processes,  solutions  of  dX  +  aX  dt  =  dL  .  It 

c 

may  be  useful  to  have  the  law  of  probability  of  the  solution  X. 


Corollary  IL.2.6:  The  characteristic  function  ij>^  of  the  solution  X  is  Vu, 


t  c  R 


=  exp  (u) 


^  z  >0 


e  ^^]dn(z)dT 


Proof:  As  uX(t)  =  N(u  Ij^^^^^e  ^),  we  use  Theorem  II. 2. 3.  q 

In  Proposition  II.  1.2,  we  had  L..  non-centered  Levy  process  with  stationary 

^  2 

increments.  The  hypothesis  /  - ~  dn(z)  <  “>  is  equivalent  to 

|zl>0  1+z 

2 

/  dn(z)  <  ®  and  /  z  dn(z)  <  ®.  The  process 
|zl>l  <\z\s^ 

l/(t)  =  lira  /  [z  -  ~^^'2  dn(z)]  is  defined  and  the  new  hypothesis 

p^  p>]z|>Jl  ^  l+z 

2 

/  z  dn(z)  <  «>  implies  that  /  z  dn(z)  <  ®  for  the  bounded  measure  on 
lzl>l  Izi^ 


Iz  r  R  ;lz|>J),  then  the  decomposition 
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L*(t)  =  lim  /  z  W  (t)-tz  dn(z)]  +  t  /  dn(z) 

I>-^P>12|>1  |z|>l  1+z^ 


gives  E(L'^(t))  =  t  / 


dn(z) 


z  >.1  1  +  z 


The  process  L‘(t)  =  lim  /  [z  W.  (t)  - ^  dn(z)]  is  of  any  order  and 

p->oo  ]>^|2|>—  1  +  z^ 

2 

E(L  (t))  =  t  J  z  dn(z)  without  any  new  assumption, 

o<rzi<i 


Theorem  1 1, 2, 6:  If  / 


z  >0  1+z 


2 

dn(z)  <  “  and  if  /  z  dn(z)  <  the  solution 


X  of  the  equation  dX  +  aX  dt  =  dL  with  a  >  0,  is  a  process  with  characteristic 

function  satisfying  Vu,  t  e  F 

,  A  iu  X(t).  ,  ,  . 

‘^X(t)  ^  )  =  exp 

“  I  -  5Tr  *  /*/,  ,  l«xp  e-=‘(*-'>ldn(z)dY 

'•  ^  -«  |zl>0  1+z 

where  y  and  6  are  two  constants. 

Proof:  By  theorem  II. 1.4  for  every  t  of  R 

L(t)  =  Y  +  6  W(t)  +  L^(t) 

2 

with  the  hypothesis  /  z  dn(z)  <  <*>  ,  the  process  L  is  second  order  and  L  (t)  = 

T>  ^ 


.  ,(t)  +  t  — — y  dn(z).  Then  the  solution  X  satisfies 

^  z|>0  1+z 

X{t)  =  ^  +  <5  e"^‘’^‘’^^dW(T)  +  /^e'®*^^'’^^dL^(T) 

-00  -00 

+  2  dn(z) 


The  four  processes  of  the  second  member  are  independent.  The  characteristic 
function  of  each  one  is  known  then 


A  random  variable  Y  has  indefinitely  divisible  law  if  its  characteristic 
function  satisfies  V  R 

2 

1  j.,-  ^  2  f  ,  iuz  ,  iu2  .  ,  r,  - 

log  <j)(u)  =  ium--:^a+J  (e  -1 - )d  )(z) 

2  |2l>0  1+z^  ^ 

where  m,a  are  constants  and  where  the  nondecreasing  function  ^  is  such  that 
lim  Jfz)  =  0,  lim  }^(z)  =  0,  J^^z^d  )|(z)  <  «>  and  /q^z^J^Cz)  <  “  (B.V.  Gnedenko 

and  A.N.  Kolmogorov  [8]).  For  every  t  of  R  the  random  variables  L(t)  is 
indefinitely  divisible. 

III.  SOLUTION  OF  THE  EQUATION  DX  =  dL  IN  STRONG  DISTRIBUTION  MEANING 

The  solution  is  known  with  forcing  term  L,  a  Gaussian  measure  (D.  de 
Brucq  and  C.  Olivier  [5]).  It  is  not  difficult  to  take  for  L  a  Poisson  measure 
defined  in  J.  Neveu  [9].  We  generalize  as  much  as  possible  introducing 
IL  -measure,  we  should  say  Levy-measure  but  the  expression  is  used  with  other 
acceptance  (Proposition  1 1. 1.2). 

We  consider  the  space  L^(r'^  .R'^.dA)  where  \  is  a  positive  measure  on 
(r”  The  Lebesgue  measure  will  be  denoted  dA  =  di.  Fourier  transform 

is  used  in  the  theory  and  the  space  L^(r'’  .R'^.dA)  is  composed  of  functions 
with  complex  values. 

Definition  III. 1.1:  Let  n  be  a  positive  measure  on  (R  ,R)  with  /z^dn(z)< 

A  IL  -measure  is  a  process 


IL  =  (n,A.P^,CX(f))  ,C  ,C) 

L^(r"  ,R",dA) 


such  as 
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(a)  f  >  L^(R^  .r’^.cIA)  •>■  IL  (f)  c  L^{.'i,A,P  )  is  linear 

3 

(b)  lij(f)  =  E(e^  IL  (f)^  _  j  j  ^  ^  f(T)-l-i  zf(z)]dn(z)dX(T) 

r"  |z|>0 

The  characteristic  function  of  the  process  IL  is  given  by  b. 

Lenuna  1 1 1  ■  1 . 2 :  The  complex  function  f  L^(R^  ,R^,dX)  ->■  i|)(f)  is  continuous. 
Proof:  We  have 

k(f)l  ^  lf(T)l^dX(T)  /  lz|^dn(z) 

R  1 z I >0 

Proposition  1 1 1. 1.3:  The  IL  -measure  is  centered  and  the  covariance  F  is  the 
function  Vf,g  e  L^(R^  ,R^,dX) 

r{f,g)  =  E(IL(f)  TLliT)  =  /  z2dn(z)/fmg(T)dX(T) 

|zl>0 

Proof:  4^(u  f)  has  derivatives  at  first  and  second  order  in  u  R  and 

tjj(u  f)  =  -  ^^  /  z^dnCz)  /  |f(T)(^dX(T)  +  o(u^) 

|zl>0  r" 

Then  E(IL  (f))  =  0  and 

E  (lE  (f)l^)  =  /  z^dn(z)  J  lf(T)l^dX(T) 

1 z 1 >0  R 

with  the  linearity  of  E  ,  we  obtain  the  covariance  F  . 

Ill, 2  Equivalent  E  -Measure  on  (S’,S') 

We  use  a  theorem  for  cylindrical  measure  (A.  Badrikian  [1]).  We  restrict 
the  E  -process  to  the  nuclear  and  countably  semi-normed  space  S,  dense 
sub-space  of  L^(R^  .R’^.dx),  here  we  consider  the  Lebesgue-measure  dX(T)  =  di. 
The  characteristic  function  =  exp  of  the  linear  process  E  is  continuous 
in  t  =  0,  Then  exists  a  probability  P  on  (S’, S')  such  as  E  is  equivalent 


to  the  process 


(S'. S'.  P,  (<*,f>)„  ,  (j:,  S) 


for  every  t  of  S  the  random  variable  is  T  c  S'  ^  <T,f>  where  <,>  is  the 


duality  between  S  and  S',  its  topological  dual.  As  IL  (f)  and  <*,f>  have  the 
same  probability  law,  Ep(<*,f>)  =  0  and  Vf,g  e  S 

=  r(f,g)  =  /  2^dn(z)/^f(T)g(T)dT 

P  1 2 1 >0  f" 

2 

The  application  f  e  S  ->•  <*,f>  e  L  (S',S',P)  is  linear  and  continuous  for  the 
topology  of  ,R'^,dT).  We  denote  /fdL  the  unique  extension  of  this  appli¬ 

cation  from  the  Hilbert  space  L^(F^  ,R’',dT)  to  the  Hilbert  space  L^(S',S',P3. 
For  Borelian  B  of  r",  the  application  1^  e  L^Cf"^  ,R",dT)  ->  Jl^dL  e  L^CS',S',P) 
is  a  vectorial  measure. 

We  suppose  now  that  the  IL  -measure  is  the  process 
IL  =  (S',S',P,  (jfdL)^,  I  C). 

Instead  of  a  general  probability  space  (12, A, P  )  the  continuity  of  <})  on  S  and 

cL  Ll,^ 

the  theorem  on  cylindrical  measures,  specify  (f2,A,P  )  into  (S',S',P). 

cl  — 

III. 3  Expression  of  the  Solution  X 

The  Fourier  transform  is  defined  in  the  introduction.  We  consider 
equation  DX  =  dT  in  distribution  meaning  V<p  e  S  <DX,<P>  =  <T,4)> 

We  limit  D  to  be  linear  operators  such  a^ 

(a)  F*D  =  AF 

(b)  A  has  derivatives  of  any  orders 

(c)  multiplication  by  A  is  a  linear  operator  on  S 

(d)  the  closed  set  F  =  {x  e  f"  ;A(x)=0}  is  F.ebesgue  almost  surely  null. 

Let  0  =  {t  €  F^  ;A(t)  /  0}  the  complement  of  F  and  let  T  =  F(P(0))  where  OCO) 
is  the  distribution  space  of  the  functions  with  derivatives  of  any  order  and 
with  compact  supports  in  0.  We  observe  that  T  is  dense  in  L^(F^  ,R”,dT) . 

We  introduce  the  definition  (D.  De  Brucq  and  C.  Olivier  [S]). 


Definition  III. 3.1:  A  solution  X  of  DX  =  dT  in  strong  distribution  meaning, 
is  a  process 
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X  =  (S', S',  P,  (X(f))^  (E,C)  such  as 
Vf  e  T  DX(f)  =  <T,f>  P  -  almost  surely. 

We  obtain  the  very  general  theorem 

Theorem  I I I. 3. 2:  For  any  probability  P  on  (S', S')  the  solution  X  in  strong 
distribution  meaning  of 

DX  =  dT 

with  D  satisfying  properties  a),  b) ,  c) ,  d)  is  linear  and  given  by 
Vf  e  T  X(f)  =  <T,  fI 


F  f  ^ 


Proof:  For  f  in  T  =  F*(P(0)),  the  function  F 

fF-f 


Ff 

-j-  1  is  also  in  T  dense  subspace 


of  S.  Then  VT  e  S'  <  T,  F 
of  D,  then 

DX(f)  =  X(D*f)  =  <T 


>  is  defined.  Moreover  we  note  D*  the  adjoint 


=  <T,F  ~  >  =  <T,TFf>  =  <T,f> 


When  the  process 

(S',S',P,  (<*,f>)j,C,C  >  is  additive 

we  note  <T,f>  =  /fdT.  With  IL -measure,  the  characteristic  function  of  the 
solution  X  is  known; 


Theorem  III. 3. 3:  The  characteristic  function  of  the  solution 
X  =  (S',S',P,  (f¥dL)j,<!:,C) 

of  the  equation  DX  =  dL  with  the  prior  hypothesis  on  the  linear  operator  I) 
and  on  the  IL  -measure,  is  given  by 


log  4^ 


X(f) 


■  r" 


z  >0 


[exp 


1  zr 


-.-iz 


]dn(z)dt 


Proof:  The  solution  X  is 
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Dcfinition  III. 3. 5:  For  the  solution  X  of  the  spatio-temporal  equation 
DX  =  dL,  the  covariance  function  at  a  given  frequency  v  =  ^  is  the  function 


s(a),x,y.z)  =  — j  /  h^dn(h)/ 


i(xc+yn+z!;) 


(2ir)^  |h|>0  R-"  |A(a.,C.n,C) 


-dcdndj; 


A  I  h  I  >0^ 

We  note  c  =  — ^ — I - = -  a  normalization  factor.  We  consider  the  linear 

(2^) 

partial  derivative  operator 

A  P  k 

D  =  I  A  (^a,  .  — .)  with  p  e  R  and  a.  .  e  I  . 
k=0  j  9t^ 

We  obtained  Fd  =  AF  with 


A  =  I  la.  .(-i)Ja)^]  . 


Properties  a),  b) ,  c),  d)  for  D  are  checked  easily.  We  observe  that  A  is 

3  2  A  2  2  2 

invariant  by  rotations  of  ,  and  is  function  of  k  =  .  Then 

E  2k  r  i  i 

A(a),x)  =  2,  (■!)  x  [W  ('i)  w  ]  is  algebraic  in  <  with  2p  complex  roots 

k=0  j 

P  2  2 

k=l,2,...,p  and  A(u),k)  =  1^  kSi 

1  (u))  =  (-l)P  J(X 
P  J  P,J 

In  that  case,  it  is  possible  to  perform  the  integrations  that  appear  in  the 
expression  of  the  function  s. 

4 

Theorem  I I I. 3. 7:  In  the  spatio-temporal  space  F  ,  the  covariance  function  s 
at  a  given  frequency  of  the  solution  X  of  the  equation  DX  =  dL  with 

P  ,  -,j 

D  =  I  A'^da  .  —  J  P  €  Fa  £  (i: 

k=o  j  at-* 


is  equal  to 
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r,  -b,  Si  la,  H  -la,  Si 

P  k  ,  k  k 

■  ^  c  V  6  1  e  e 

~  I  i2  ^  b  ?  '  -T 

2  1  k=0  °k  2a,  i  d,  cT, 

'  p  k  I  k  k 

u  ^  o2  2  2  2 

where  aj  5,  =  x  +y  +z 


bl  =  a^-ibj^  aj^>  0  b^  >  0  and  1  AOo.k)  1 2=  [  Ip  |  ^ 

dk  =  n  if  P  1  and  d^  =  1. 

i=i  ■’  J 

We  have  supposed  the  roots  Xj^  in  the  complex  plane  with  bj^  strictly  positive; 
it  is  the  analog  of  a  >  0  in  Proposition  11.2.4. 

Proof:  The  theorem  is  known  for  p=l  (D.  de  Brucq  (4]).  We  have  to  compute 

^iCxC+yn+zrJ 

s(u),x,y,z)  =  c  /  - 5 - dr.dnd!; 

R  1A(u),k)P 

2  1 2  P  2  2  -2 

with  lA(fjj,K)|  =  ll  I  fl  (k  -X,)(<  -Xn.  We  used  spherical  coordinates  and 
p  k=i  " 

(I, I.  Gihman  and  A.V.  Skorohod  [7]) 


,  .  1  f+«>  iKP  K 

■Oo,x,y,z)  =  c  ^  e  - - 


1a(w,k)  1' 


Let  f(K)  = 


and  J  =  J  „  f(K)d< 


P  ,  2  ,  ,2,  2  r2. 


we  perform  the  integration  using  residual  method.  The  p  poles  are  strictly 
complex  by  hypothesis  and  =  a^^+ibj^,  -Aj^  =  -aj^+ibj^  k=l,2,...,p  are  the 
poles  in  the  upper  half  plane.  Then 


s  =  2iR  ^  [Res(f )+Res(f ,-A,  )] .  We  have 
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iX  I 

Res(f.X|^)  =  - 5 - 

lal  -  xJ)(a2  -  xlnil 


as 


n 

?=i 

2,j«k 


Sia^bk 


A,  =  4ia  b, 
k  k  k 


e  1 


‘‘k  *“k\ 


and  also 


-iA,  i 

k  , 

Res(f.  -AJ=  i - 

d,  -8ia,  b, 
k  k  k 


Going  back  to  the  expression  of  s,  we  find 

.  ^  1  1 
s(a),x,y,z)  =  c  -  - j  J 

nil  li  1^ 

'  P 


-L-2inI  ~± 


ij  j2  niH  k=l  Sia^bk 

P 


.  l1 


For  Gaussian  measures,  for  Poisson  measure,  results  of  Ill.l  and  III. 2  are 
valid  and  we  have  equivalent  processes  on  (S',S',P).  These  processes  are 
linear  and  for  every  f  in  .R^^jdr)  the  characteristic  functions  are 

given  by 

log  =  -I  J  |fl^dT  and  log  =  /  (exp  if-l-if)di 

F  F 
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Thcsc  centered  processes  and  the  FL -measures  have  the  same  correlation  function, 
hvery  forcing  terms  of  these  types,  give  solution  X  with  the  same  second  order 
properties  1  Moments  of  greater  order  are  necessary  to  make  differences  between 
the  forcing  terms,  for  p  =  1 


s(a),  ?;) 


sin  ki! 


where  =  k-iy  is  the  complex  wave 


number.  Direct  verification  in  deep  se.i  water  of  this  formula  is  factible: 

k((n)  is  the  wave  number  at  pulsation  'o, 

,  ((!))  is  the  dumping  term  at  pulsation  w. 

For  two  points  at  a  distance  9.,  the  correlation  of  the  filtered  observations 
at  frequency  v  of  the  pressure  p  is  the  function  s(ai,£).  The  main  assumption 
to  obtain  the  result  is  that  p  satisfies  any  equation 

with  I,  Gaussian,  Poisson,  or  any  centered  Levy-measure! 
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